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Resource Allocation and Rate Gains 
in Practical Full-Duplex Systems 

Jelena Marasevic, Jin Zhou, Harish Krishnaswamy, Yuan Zhong, Gil Zussman 


Abstract —Full-duplex communication has the potential to sub¬ 
stantially increase the throughput in wireless networks. However, 
the benefits of full-duplex are still not well understood. In this 
paper, we characterize the full-duplex rate gains in both single¬ 
channel and multi-channel use cases. For the single-channel 
case, we quantify the rate gain as a function of the remaining 
self-interference and SNR values. We also provide a sufficient 
condition under which the sum of uplink and downlink rates on 
a full-duplex channel is concave in the transmission power levels. 
Building on these results, we consider the multi-channel case. For 
that case, we introduce a new realistic model of a compact (e.g., 
smartphone) full-duplex receiver and demonstrate its accuracy 
via measurements. We study the problem of jointly allocating 
power levels to different channels and selecting the frequency 
of maximum self-interference suppression, where the objective is 
maximizing the sum of the rates over uplink and downlink OFDM 
channels. We develop a polynomial time algorithm which is nearly 
optimal in practice under very mild restrictions. To reduce the 
running time, we develop an efficient nearly-optimal algorithm 
under the high SINK approximation. Finally, we demonstrate 
via numerical evaluations the capacity gains in the different use 
cases and obtain insights into the impact of the remaining self¬ 
interference and wireless channel states on the performance. 

Index Terms —Full-duplex, modeling, resource allocation. 

1. Introduction 

Full-duplex (FD) communication - simultaneous transmis¬ 
sion and reception on the same frequency channel - holds great 
promise of substantially improving the throughput in wireless 
networks. The main challenge hindering the implementation 
of practical FD devices is high self-interference (SI) caused 
by signal leakage from the transmitter into the receiver circuit. 
The SI signal is usually many orders of magnitude higher than 
the desired signal at the receiver’s input, requiring over lOOdB 
(i.e., by 10^^ times) of self-interference cancellation (SIC). 

Cancelling SI is a very challenging problem. Even though 
different techniques of SIC were proposed over a decade ago, 
only recently receiver designs that provide sufficient SIC to 
be employed in Wi-Fi and cellular networks emerged (see 
Q and references therein for an overview). Exciting progress 
was made in the last few years by various research groups 
demonstrating that a combination of SIC techniques employed 
in both analog and digital domains can provide sufficient SIC 
to support practical applications 0 (T§ 

While there has been significant interest in ED from both 
industry and academia p|-p4|, the exact rate gains resulting 
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Eigure 1: Some possible uses of full-duplex: (a) simultaneous 


UL and DL for one MS;[^UL and DL used by two different 
MSs and caused inter-node interference (red dashed line), 
simultaneous UL and DL over OEDM channels. 


from the use of FD are still not well understood. The first 
implementations of ED receivers optimistically envisioned 
100% rate improvement (e.g., j^, fT^). To achieve such an 
increase in data rates, the ED receiver would need perfect SIC, 
namely, to cancel SI to at least one order of magnitude below 
the noise floor to render it negligible. The highest reported 
SIC |T^ , however, suppresses the SI to the level of noise. 

Despite this insufficient cancelling capabilities, much of the 
work on ED rate improvement assumes perfect SIC in the 
ED receiver |T^-pT| . While non-negligible SI has also been 
considered ||22|-||24|, there are still no explicit bounds on the 
rate gains for given FD circuit parameters and parameters of 
the wireless signal. Moreover, from a modeling perspective, 
the frequency selectivity of SIC has not been considered in 
any analytical work. This is an important feature that is 
inherent in conventional compact implementations of an FD 
receiver, such as that found in small-form factor mobile devices 
(e.g., smartphones and tablets), where frequency selectivity is 
mainly a consequence of the cancellation in the RE domain Q 

A. Summary of Contributions 

The main contribution of this paper is a thorough analytical 
study of rate gains from FD under non-negligible SI. We 
consider both single-channel and multi-channel orthogonal 
frequency division multiplexing (OFDM) scenarios. Eor the 
multi-channel case, we develop a new model for frequency- 
selective SIC in small-form factor receivers. Our results pro¬ 
vide explicit guarantees on the rate gains of ED, as a function 
of receivers’ signal-to-noise ratios (SNRs) and SIC profile. Our 
analysis provides several insights into the structure of the sum 
of uplink (UL) and downlink (DL) rates under ED, which will 
be useful for future work on ED MAC layer algorithm design. 

Specifically, we consider three different use cases of ED, 
as illustrated in Eig. (i) a single channel bidirectional link, 

^See our recent work |25) , and Section |n^ for more details. 
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where one mobile station (MS) communicates with the base 
station (BS) both on the UL and on the DL (Fig. |l|[a)|); 


(ii) two single channel unidirectional links, where one MS 
communicates with the BS on the UL, while another MS 
communicates with the BS on the DL (Fig. |l||^b)|); and (hi) a 


multi-channel bidirectional link, where one MS communicates 
with the BS over multiple OFDM channels, both on the UL 
and on the DL (Fig. |l|[c)| ). 

1) Models of Residual SI: For SI, we consider two 
different models. For the BS in all use cases and the MS in 
use case (i), we model the remaining SI after cancellation as 
a constant fraction of the transmitted signal. Such design is 
possible for devices that do not require a very small form 
factor (e.g., base stations), and was demonstrated in 

In the multi-channel case, we rely on the characteristics 
of RFIC receivers that we recently designed | [25| , | [26| and 
develop a frequency selective model for the remaining SI in 
a small form-factor device (Section We demonstrate the 
accuracy of the developed model via measurements with our 
receivers p5| , | [26| . We note that a frequency-selective profile 
of SIC that we model is inherent to RF cancellers with flat 
amplitude and phase response (see Section |IIJ. A mixed-signal 
SIC architecture GD where the digital TX signal is processed 
and upconverted to RF for cancellation does not necessarily 
have flat amplitude and phase response. However, we do not 
consider this architecture because it requires an additional up- 
conversion path compared to the architecture of this work, and 
this additional path introduces its own noise and distortion, 
limiting the resultant RF SIC. 

2) Sum Rate Maximization: We focus on the problem of 
maximizing the sum of UL and DL rates under FD (referred 
to as the sum rate in the rest of the paper). This problem, 
in general, is neither concave nor convex in the transmission 
power levels, since the remaining SI after cancellation depends 
on the transmission power level. Due to the lack of a good 
problem structure, existing analytical results (see e.g., 
p4| ) are often restricted to specialized settings. Yet, we obtain 
several analytical results on the FD rate gains, often under mild 
restrictions, by examining closely the structural properties of 
the sum rate function. 

Single-Channel Results. In the single-channel cases, we prove 
that if any rate gain can be achieved from FD, then the gain 
is maximized by setting the transmission power levels to their 
respective maximum values. This result is somewhat surprising 
because of the lack of good structural properties of the sum 
rate. We then derive a sufficient condition under which the 
sum rate is biconcav^in both transmission power levels, and 
show that when this condition is not satisfied, one cannot 
gain more than Ib/s/Hz (additively) from FD as compared to 
time-division duplex (TDD). We note that although the model 
for the remaining SI in the single channel case is relatively 
simple, it nonetheless captures the main characteristics of the 
FD receivers. Moreover, the results for the single channel case 
under this model are fundamental for analyzing the multi¬ 
channel setting, and often extend to this more general setting. 


function is biconcave, if there exists a partition of variables into two 
sets, such that the function is concave when variables from either set are fixed. 


Multi-Channel Results. In the multi-channel case, we use 
the frequency-selective SI model for the MS receiver that is 
introduced in Section IV-A and motivated by FD implementa¬ 
tion challenges discussed in Section |I^ Based on this model, 
we study the problem of transmission power allocation over 
OFDM channels and frequency selection, where the objective 
is to maximize the sum of the rates over UL and DL OFDM 
channels (in this case, frequency refers to the frequency of 
maximum SIC of the SI canceller). Although in general it is 
hard to find an optimal solution to this problem, we develop 
an algorithm that converges to a stationary point (in practice, 
a global maximum) under two mild technical conditions. 
One condition ensures that the sum rate is biconcave in 
transmission power levels. This restriction is mild, since we 
prove that when it does not hold, the possible gains from 
FD are small. The other condition imposes bounds on the 
magnitude of the first derivative of the sum rate in terms of 
maximum SIC frequency, and has a negligible impact on the 
sum rate in OFDM systems with a large number of channels, 
because it can only affect up to 2 OFDM channels (see Section 
|VI-A| for more details). 

Although the algorithm in practice converges to a near- 
optimal solution and runs in polynomial time, its running 
time is relatively high. Therefore, we consider a high SINR 
approximation of the sum rate, and derive fixed optimal power 
allocation and maximum SIC frequency setting that maximizes 
the sum rate up to an additive e in time 0(iT log (1/e)), for 
any given e, where K is the number of channels. 

Numerical Results. Finally, we note that throughout the paper, 
we provide numerical results that quantify the rate gains 
in various use cases and illustrate the impact of different 
parameters on these gains. For example, for the multi-channel 
case, we evaluate the rate gains using measured SI of our RFIC 
receiver p5| , p6| . We use algorithms for the general SINR 
regime and for the high SINR regime and compare their results 
to those obtained by allocating power levels equally among 
the OFDM channels. Our results suggest that whenever the 
rate gains from FD are non-negligible, all considered power 
allocation policies yield similar rate gains. Therefore, one of 
the main messages of our work is that whenever it is beneficial 
to use FD, simple power allocation policies are near-optimal. 


B. Organization of the Paper 

The rest of the paper is organized as follows. Section 
reviews related work and Section |m] outlines the challenges 
in implementing FD receivers. Section [Tvl introduces the new 
model of a small form factor FD receiver, and the model for the 
various use cases. Section [V| provides analysis and numerical 
evaluation for the sum rate maximization on a single channel 
for use cases (i) and (ii). Sections [Vl| and [Vll| provide analysis, 
algorithms, and numerical evaluation for use case (iii). We 
conclude in Section IVIIII 

II. Related Work 

Possible rate gains from FD have been studied in 
p4| , with much of the work p~8|- (D focusing on perfect 
SIC. Unlike this body of work, we focus on rate gains from 
FD communication under imperfect SIC. 
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Figure 2: Block diagram of a full-duplex transceiver employ¬ 
ing RF and digital cancellation. 

Non-negligible SI has been considered in p^-p4|. A 
sufficient condition for achieving positive rate gains from FD 
on a bidirectional link has been provided in p2| , for the 
special case of equal SINRs on the UL and DL. This condition 
does not quantify the rate gains. 

Power allocation over orthogonal bidirectional links was 
considered in and for MIMO and OFDM systems, 
respectively. The model used in p4| assumes the same amount 
of SIC and equal power allocation on all channels, which is a 
less general model than the one that we consider. 

A more detailed model with different SIC over OFDM 
channels was considered in . The model from does not 
consider dependence of SIC in terms of canceller frequency 
(although, unlike our work, it takes into account the transmit¬ 
ter’s phase noise). Optimal power allocation that maximizes 
one of the rates when the other is fixed is derived for equal 
power levels across channels, while for the general case of 
unequal power levels, | [23] only provides a heuristic solution. 

Our work relies on structural properties of the sum rate 
to derive near-optimal power allocation and maximum SIC 
frequency setting that maximizes the sum rate. While the 
model we consider is different than p2| , | [2^ , we provide 
a more specific characterization of achievable rate gains, and 
derive results that provide insights into the rate dependence on 
the power allocation. These results allow us to solve a very 
general problem of rate maximization. 


III. FD Implementation Challenges 


In this section, we overview the challenges associated with 
the implementation of compact FD radios. These challenges 
motivate the model of remaining SI that is introduced in 
Section [TV- A| and used in the design of sum-rate maximization 
algorithms (Section VI ). 

Fig. shows the block diagram of a full-duplex transceiver. 
There are two antenna interfaces that are typically considered 
for full-duplex operation: (i) an antenna pair and (ii) a circu¬ 
lator. The advantage of using a circulator is that it allows a 
single antenna to be shared between the transmitter (TX) and 
the receiver (RX). SIC must be performed in both the RF and 
digital domains to achieve in excess of lOOdB SI suppression. 
The RF canceller taps a reference signal at the output of the 
power amplifier (PA) and performs SIC at the input of the 
low-noise amplifier (LNA) at the RX side ( 27 ) 

Typically, 20-30dB of SIC is required from the RF, given 
that the antenna interface typically has a TX/RX isolation of 
20-30dB |[28l. Thus, an overall 50-60dB RF TX/RX isolation 





(b) 


Figure 3:|(a)|RFIC receiver with RF SI cancellation |25| , 
and the two antenna interfaces used in our measurements: 
an antenna pair and [(c)| a circulator. 


is achieved before digital SIC is engaged. This amount of 
RF TX/RX isolation is critical to alleviate the RX linearity 
and the analog-to-digital conversion (ADC) dynamic range 
requirements 0, p7| . Digital cancellation further cancels the 
linear SI as well as the non-linear distortion products generated 
by the RX or the RF canceller. 

A mixed-signal SIC architecture has been proposed in (T§, 
where the digital TX signal is processed and upconverted to 
RF for cancellation. However, this requires a separate up- 
conversion path which introduces its own noise and distortion. 
Moreover, the noise and distortion of the TX analog and 
RF circuits (such as the power amplifier) are not readily 
captured in the cancellation signal, limiting the resultant RF 
SIC. In addition, the dedicated up-conversion path results in 
area and power overhead. Because of these reasons, we are 
not considering this SIC architecture in this paper. 

For wideband SIC, the transfer function of the canceller 
must closely track that of the antenna interface across fre¬ 
quency. However, the frequency dependence of the inherent 
antenna interface isolation together with selective multi-path- 
ridden SI channels render this challenging for the RF canceller 
in particular. The net antenna interface isolation amplitude 
and phase response can vary significantly with frequency. A 
rapidly-varying phase response is representative of a large 
group delay, requiring bulky delay lines to replicate the 
selectivity in the RF canceller |T^ , p7| . 

The fundamental challenge associated with wideband SIC 
at RF in a small form-factor and/or using integrated circuits 
is the generation of large time delays. The value of true time 
delay is linearly proportional to the dimension of the delay 
structure and inversely proportional to the wave velocity in 
the medium. To generate Ins delay in a silicon integrated 
circuit, a transmission line of 15cm length is required as the 
relative dielectric constant of silicon oxide is 4. A conventional 
integrated RF SI canceller with dimensions less than Imm^ 
will therefore exhibit negligible delay. Note that the canceller 
phase response can be calculated by integrating the delay 
with respect to frequency, and conventional integrated RF 
SI cancellers typically have a fiat amplitude response flSl , 
p6| . Therefore, the amplitude and phase response of the 
canceller can be assumed to be flat with respect to frequency 
when compared with antenna interface isolation, limiting the 
cancellation bandwidth 0 , (26) 

While achieving wideband RF SI cancellation using in¬ 
novative RFIC techniques is an active research topic (e.g., 
frequency domain equalization based RF SI cancellation in 
p9|), in this paper we focus on compact fiat amplitude- and 
phase-based RF cancellers, such as the one we implemented 
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Figure 4: Measured isolation amplitude and group delay of|(a) 
a PCB antenna pair and |(b)| a commercial 2110-2170 MHz 
miniature circulator from Sky works p8| , and the resultant 
TX/RX isolation using the integrated RF canceller with flat 
amplitude and phase response from |[25|, p6| with the 


antenna pair and (d) the circulator compared to the SIC model. 


in the RFIC depicted in Fig. |^|^a)| | [25| , | [26| . 

In 1251 and l26| , the RF canceller is embedded in the RX’s 
LNA, and consists of a variable amplifler and a phase shifter. 
The RF canceller adjusts the amplitude and the phase of a 
TX reference signal tapped from the PA’s output performing 
SIC at the RX input. Thanks to the co-design of RF canceller 
and RX in a noise-cancelling architecture, the work in |25l 
and l26l is able to support antenna interface with about 20dB 
TX/RX isolation with minimum RX sensitivity degradation. 

We measured isolation amplitude and group delay response 
of (i) a PCB antenna pair (see Fig. and (ii) a commercial 
2110-2170MHz miniature circulator from Sky works m (see 


Fig. |^[c)| ). The results are shown in Fig. |^|^a)| and Fig. HE 


respectively. The resultant TX/RX isolations using an RF 
canceller with flat amplitude and phase response after the 
antenna interfaces (i) and (ii) are shown in Fig. Nfr)| and 
Fig. |4[d)[ respectively. As Fig. |4[c)| and Fig. |4[d)| suggest, for 
-60dB TX/RX isolation after RF cancellation, the bandwidths 
are about 4MHz and 2.5MHz, respectively. 


IV. Model 

We consider three use cases of FD: (i) a bidirectional link, 
where one mobile station (MS) communicates with the base 
station (BS) both on the UL and on the DL (Fig. |l|[a)| ), (ii) 
two unidirectional links, where one MS is communicating with 
the BS on the UL, while another MS is communicating with 
the BS on the DL (Fi g. |l|[b)| ), and (iii) multiple orthogonal 
bidirectional links (Fig. [JcJ^ Note that in (ii) only the BS is 
operating in FD. 

For the multi-channel FD (use case (hi)), we assume that the 
network bandwidth of size B is subdivided into K orthogonal 
frequency channels of width B/K each, and index the fre¬ 
quency channels with k G iT}. An example of such 

sub-channelization is OFDM with each frequency channel 
consisting of an integral number of subcarriers. 


For all notation that relates to the BS, we use b in the 
subscript. For the notation that relates to the MS in use cases 

(i) and (iii), we use m in the subscript, while in the use case 

(ii) we use mi and m 2 to refer to MS 1 and MS 2, respectively. 
Summary of the main notation is provided in Table 

The transmission power of a station u G {6, m, mi, m 2 } on 
channel k is denoted by Pu,k^ where k G AT}. In use 

cases (i) and (ii), k is omitted from the subscript, since we 
consider a single channel. 


A. Remaining SI 


Single-channel FD. For single-channel FD, we assume that 
the remaining SI both at the BS and at an MS can be expressed 
as a constant fraction of the transmitted power. In particular, 
if the BS transmits at the power level P^, the remaining SI 
is RSIh = QbPb, where gb is a constant determined by the 
hardware. Similarly, if an MS transmits at the power level 
Pm, its remaining SI is RSIm = QmPm• 

Multi-channel FD. We assume that the FD receiver at the BS 
has frequency-flat SIC profile, meaning that the remaining SI 
at the BS on channel k is RSIb^k = dbPb.k, where gb is a 
constant. We note that such FD receiver design is possible to 
implement in devices that do not require small form factor of 
the circuit (e.g., a BS or an access point (AP)), and has been 
reported in (D 

In the rest of this section, we describe the mathematical 
model of the remaining SI for a small form factor device 
(MS). We consider a compact/RFIC FD receiver with a 
circulator at the antenna interface, described in Section ||n| 
and assume a frequency-flat amplitude and phase response of 
the canceller, denoted by \Hc^r\ and ZHc,r, respectively. The 
amplitude and phase response of the canceller are assumed to 
be programmable but constant with frequency. 

For the antenna interface’s TX/RX isolation, we assume a 
flat amplitude response \HA{f)\ = const = \Ha\ and a con¬ 
stant group delay equal to r, so that HA{f) = 

(recall that the measured amplitude and group delay response 
are shown in Fig. For the digital SIC, denoted hy SICd, 
we assume that the amount of cancellation is constant across 
frequency, as delay can be easily generated in the digital 
domain. Let fk denote the central frequency of the channel, 
so that fk = fi + {k — PjBjK. Then, the remaining SI after 
cancellation can be written as: 


=\PmA^A - Hc,r)SIC-^^\ 

=P„,fe|(|Ff^|2 + \Hc,r\^ - 2\Ha\\Hc,r\ 

■cos{ZHA{fk) + ^Hc,R))\SIC^\ ( 1 ) 

Note that in Pm,k is the MS transmission power on 
channel k, Pm,k{HA — Hc,r) is the remaining SI after the 
RF SIC, and Pm^ki^A — Hc,R)SICf)^ is the remaining SI 
after both the RF and digital SIC. 

We assume a common oscillator for the TX and RX, with 
the phase noise of the oscillator being good enough so that it 
does not affect the remaining SI. 

The RF canceller’s settings can be programmed in the held 
to adjust the frequency at which peak SIC is achieved flSl , 
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Table I: Nomenclature. 


m 

Subscript notation for an MS 

b 

Subscript notation for the BS 

K 

Total number of OFDM channels 

k 

Channel index, k G {1, ..., K} 

U, V 

Station indexes, u,u G {5,m,mi,m 2 } 

Pu,k 

Transmission power of station u on channel k 

Pu 

Maximum total power: Pn,/e < Pu 

9h 

Remaining SI at the BS per unit transmitted power 

c 

Position of the maximum SIC frequency, c G M 
Remaining SI at an MS: (i) per unit transmitted power 

Qm 

huv,k 

for P = 1; (ii) per unit transmitted power at unit 
distance from the maximum SIC frequency for P > 1 
Wireless channel gain for signal from u to u on channel 

k, for u ^ u 

Nu 

Thermal noise at station u 

^uv,k 

SNR of signal from u to u on channel k, where u ^ v 

^uu,k 

XINR at station u, channel k 

max 

luv,k 

yuVjk for Pu — Pu 

Vk 

Total rate on channel k 

r 

Sum rate: r = 


With the amplitude {\Hc^r\) and the phase (ZHc^r) of 
the RF canceller set to \Ha\ and —ZHA{fc)^ respectively, 
peak SIC is achieved at frequency /c. Therefore, the total 
remaining SI at the MS on channel k can be written as: 


RSI^^k = 2\HAfPm,kil - cos(27rT(/fe - fc)))SIC 




where r is the group delay from the antenna interface with 
a typical value at the order of Ins (which agrees with the 
measured group delay in Fig. Frequency bands used by 
commercial wireless systems are at most 10s of MHz wide. 
It follows that 27rr(//c — fc) « 1, and using the standard 
approximation cos(x) « l — x^/2 for x « 1, we further get: 

RSIm,k ~ \HA?PmA‘^TrTfUk - fc?SIC-^\ 

Recalling that //c = /i + — 1)B/K = /o + kB/K for 

fo = fi - B/K, and writing fc as /^ = /o + cB/K, for 
c G M, we can combine all the constant terms and represent 
the remaining SI as: 

RSIm,k — QmPm^kiJ^ P) 5 (^) 


where Qm = \Ha\^{‘ 2.11 r^^B/SIC. Note that even 
though in this notation we allow c to take negative values, 
we will later show that in any solution that maximizes the 


sum rate it must be c G (1, Ff) (Lemma 6.3). 

Fig. ^I^d) shows the TX/RX isolation based on Eq. ^ and 
based on measurement results. The parameter gm in Eq. 0 
was determined via a least square estimation. The modeled 
TX/RX isolation based on Eq. Q is also compared to the 
measured TX/RX isolation of the canceller with the antenna 
pair interface in Fig. ^|^c) As Fig. shows, our model of 
the remaining SI closely matches the remaining SI that we 
measured with the RFIC FD receiver presented in |[26|. 


B. Sum Rate 

The total transmitted power of each station is assumed to 
be bounded as follows. In use cases (i) and (ii): P 5 < Pb, and 
each Pm, Pmi , Pm 2 ^ Pm- In use case (hi): Pu,k < Pu, 

where u G Pu > 0. The channel gain from station u 


to station v on channel k is denoted by huv,k in use case (hi) 
and by huv in use cases (i) and (ii). The noise level at station 
u is assumed to be equal over channels and is denoted by Nu. 
We assume that the channel states and noise levels are known. 

For the signal transmitted from u to v, where u^v G 
{ 6 ,m,mi,m 2 }, u V, and either u = b or v = b, we 
let juv,k = denote signal to noise ratio (SNR) at 

V on channel k. Similarly as before, in use cases (i) and 
(ii), index k is omitted from the notation. In the use case 
(ii), 7 mim 2 denotes the (inter-node-)interference to noise ratio 
(INR). Self-interference to noise ratio (XINR) at the BS is 
denoted by 755 = in use cases (i) and (ii), and by 

766 ,/c = {^^^^- ^INR at the MS is denoted 

by Imm = 1mm,k = 

and (iii), respectively. 

We use Shannon’s capacity formula for spectral efficiency, 
and let log(.) denote the base 2 logarithm, ln(.) denote the 
natural logarithm. We use the terms “spectral efficiency” and 
“rate” interchangeably, as the spectral efficiency on a channel 
is the rate on that channel normalized by B/K. 

In use case (i), the sum rate on the channel is given as: 

r = log (1 + + log (1 + • (3) 

1 + 766 2 V 1 -h 7mm 2 


Observe that 17 ^ and are signal to interference-plus- 

noise ratios (SINRs) on the UL and DL, respectively. We will 
refer to Vm = log (^1 + i+m) ~ 

log (1 -h I the DL rate. 

Similarly as for (i), the sum rate for use case (ii) is: 


r = 


log (^1 


7mi6 

1 + 766 


) + log 


T 6 m 2 


1 P 7mim2 


)■ 


(4) 


Finally, in use case (iii), the rate on channel k is given as: 


Tk = 


= log ( 


1 + 


7 m 6 ,/c 

1 + 766,/c 


) + log ( 


1 + 


76m,/c 
1 H“ ^mm,k 


)■ 


(5) 


while the sum rate (on all channels) is r = Ylk=i 

The objective in all problems considered is to maximize 
r subject to the upper bound on total transmitted power 
and non-negativity constraints. In use cases (i) and (ii), 
the variables are P 5 and P^, while in the use case (iii), the 
variables are c, P^^, and Pm,/c, for k G {1, ...,Pr}. 

For the purpose of comparison to TDD systems, we will 
sometimes also consider TDD rates. We denote by ^ = 
log(l + 7 ™^"^) and r^^ b = 1®S(I +7^"^) maximum UL 
and DL TDD rates, respectively, where 7 ™^^ = ymb{Pm)^ 
"ibyrT ~ lhm{Pb)- The maximum achievable TDD rate can 
then be written as = maxjr^g^ m’^TDD 61 - 


V. Single Channel FD 
A. A Bidirectional FD Link 

In this section, we derive general properties of the sum rate 
function for use case (i) (Fig. m- 

First, we show that if it is possible for the FD sum rate 
to exceed the maximum TDD rate, it is always optimal for 
the MS and the BS to transmit at their maximum respective 
power levels (Lemma [5.1| ). This result is somewhat surprising, 
because in general, the FD sum rate function does not have 
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good structural properties, i.e., it need not be convex or 
concave in the transmission power variables. 

Building upon this insight, we quantify the FD rate gains 
by comparing the FD sum rate to corresponding TDD rates 
(Section |V-A2| ). More specifically, we define a metric that 
characterizes by how much the FD capacity region extends the 
corresponding TDD capacity region, and provide a sufficient 
condition on the system parameters for rate gains to hold. 

Finally, we establish a sufficient condition for the FD sum 
rate function to be biconcave in transmission power levels 


r = 


= log ( 


Nb + QbPb ) 


+ 


log ( 


1 + 


9mPrr 


= log (1 + 7) , where 
hfnbPm , ^bmPb 


)■(' 


hbmPb 

QmPri 


:)) 


7 = 


+ 


+ 


hmbPm 


hbmPb 


P'b + 9bPb Nm + druPm P'b + QbPb Nm + druPy 


= - 2 - 


{Nm P gmPmY \Nb P QbPb 

^bmPbdm 


■ 9m 


hmbP^rr 


dpj 

From ([^ and Q: 


+ 9mPmY ^^bP 9bPb 


■ 9m^ • 


( 6 ) 

(7) 


1) If 


> 0, i.e.. 


Nr+gbPb dP^'^ — ^ dPm 

7 is concave and strictly increasing in Pm when Pb is 

fixed, and therefore maximized for Pm = Pm- 


2) If 


hmb ^ 

Nb+gbPb 

strictly convex in Pj 


< 0, then 




> 0, i.e., 7 is 




r^/r; 


max 
TDD,6 


(Section V-A3). This condition imposes very mild restrictions 
on the XINRs at the BS and the MS. Moreover, the established 
condition extends to the multi-channel scenario (use case (iii)), 
where it plays a crucial role in deriving an algorithm for the 
sum rate maximization that converges to a stationary point 
that is a global maximum in practice (Section [VI-Bl[ ). Without 
such a condition, the problem would not have enough structure 
to be amenable to efficient optimization methods. 

1) Power Allocation: 

Lemma 5.1: If there exists an FD sum rate r that is higher 
than the maximum TDD rate, then r is maximized for Pm = 

Pm, Pb = Pb- 

Proof: From ([^, the sum rate can be written as: 

^ ^ hmbPm ^ I 1_ f I hbmPb 


(b) 

Figure 5: TDD and FD capacity regions, and FD extension. 
The capacity region is plotted for equal maximum SNRs: 
7m6^ = 7^^ = l^^mb cases of maximum XINRs: 


(a) 75 ^"" = 1 , = 1 and (b) 75 ^"" = 1 , 7 “^" = 10 . 


capacity region to the corresponding TDD capacity region. 
Let rb = log(l + 77r)’ = 1°S(1 + T+^) denote DL 


l+ 7 bf> ^ 


Taking partial derivatives of r directly does not provide con¬ 
clusive information about the optimal power levels. Instead, 
we write r as an increasing function of another function that 
is easier to analyze. Specifically: 

^mbPrr. 


and UL rates, respectively and let ^ = log(l + 5^^), 
^TDD m = lc)g(l + 7m6^) ^^^ote the maximum TDD rates. 
The FD capacity region is the set of all points (rb^Vm) such 
that Pm ^ [0,Pm], Pb ^ [0,A], while the TDD capacity 
region is the convex hull of the points (0,0), (r^g 0), 

max ’ 

and (0,r^'g,„). We also let Sb = log(l + and 

Sm = log(l + ^^^max ) be the FD DL and UL rates when 

' 'bb _ _ 

both stations transmit at their maximum power levels Pb^ Pm- 
Fig. shows FD and TDD capacity regions for symmetric 
maximum SNRs "fmt^ = 

^max _ ^max 
Jbb fmm 

the axes are normalized by ^ and respectively. To 

determine the points at the boundary of the FD capacity region, 
57] as follows. For Vb = asb, where a G 

is maximized 


7 ^^ and two cases of maximum 
XINRs: = 7 “S? = 1 and 7 ^^^^ = 1, 7 ^^;^ = 10. Here, 


we apply Lemma 
(0,1), Lemma 5.1 
for 


implies that the UL rate 


Since r is strictly increasing in 7, to maximize r it suffices 
to determine Pm^Pb that maximize 7. The first and the second 
partial derivative of 7 with respect to Pm are: 

57 _ h mb I ^bmPb f hmb^rr 

dPm Nb + 9bPb 

9^7 


Pm = Pm, regardless of the value of P 5 . Therefore, the 
DL rate is lowered from Sb to Vb = asb by lowering Pb. The 
point {vb^Vm) at the boundary of the FD capacity region is 
then determined by solving Vb = asb for Pb, and setting Vm = 
rm{Pb^ Pm)- An analogous procedure is carried out for Vm = 
aSm, where a G (0,1). We remark that FD capacity regions 
are not necessarily convex (e.g.. Fig. 

anH ^max 

%mlmb 


b) 


for 


Pmjmb 


= OdB 


OP 2 

ur^ m 

when Pb is fixed. Therefore, 7 is 
maximized at either Pm = 0 or Pm = Pm- Note that if 
Pm = 0, there is no signal on UL, in which case FD rate 
equals the maximum TDD UL rate. 

A similar results follows for Pb by taking the first and the 
second partial derivative of 7 with respect to P^. ■ 

2) Mapping Gain over SINK Regions: In this sec¬ 
tion we quantify the FD rate gains by comparing the FD 


Lemma 15.11 states that the maximizer of the FD sum rate 
is either (r^D ft, 0 ), ( 0 ,r^“„j) or (s 6 ,Sm)- In particular, to 
see whether FD operation increases the sum rate, it suffices 
to check whether Sb P Sm > 6’^tdd ml- "^bis 

motivates us to focus on the pair ( 55 , Sm) when considering by 
how much the FD capacity region extends the corresponding 
TDD capacity region. We introduce the following definition 
(see Fig. |^|^b) for a geometric interpretation). 

Definition 5.2: FD extends the corresponding TDD capacity 
region by p - 100 % if p > 0 is the smallest number for which 
is inside the TDD capacity region. 

The following lemma provides a necessary and sufficient 
condition for the capacity region extension of p • 100 %. 

Lemma 5.3: FD extends the TDD capacity region by 
p • 100 %, where p > 0 , if and only if: 


log 


1 + 




) / max \ 

log (1 + 1+7-j 


log(l + log(l + 7mr) 


= l+p. (8) 
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7“r[dB] 


7rr[dB] 


(a) 


(b) 


Figure 6: TDD capacity region extension due to FD as a 
function of SNRs for [(a)| 7 ^^"" = 1, 

irr = h = 10 . 


7 ”“ = 1 and (b) 


Proof Sketch: The proof is based on the fact that since 
p is the smallest number for which the point ( 3 ^, 
is in the TDD capacity region, must lie on the 

line connecting ^ and ^ (Hg. |^b)l), and therefore: 


j^max 

1+p ~ 'TDD,m 1+p’ 

Fig. shows the TDD capacity region extension due to 
FD operation, as a function of the received signals’ SNR, for 
BS FD receiver that cancels SI to the noise level and MS FD 
receiver that cancels SI to (i) the noise level (Fig. |^[a)| ) and (ii) 
one order of magnitude above noise (Fig. Recall from 

Definition |5.2| that the capacity region extension is computed 
for Pm = Pm and Pi, = Pi,, and therefore the differences in 
the SNRs are due to signal propagation and not due to reduced 
transmission power levels. Fig. [^suggests that to achieve non- 
negligible capacity region extension, SNRs at the MS and at 
the BS must be sufficiently high - at least as high as to bring 
the resulting SINR to the level above OdB. 

3) Sum Rate Biconcavity: In this section, we establish a 
sufficient condition for the sum rate to be (strictly) biconcave 
and increasing in Pm and in Pi, (Condition |5.4| ). We also show 
that when the condition does not hold, using FD does not 
provide appreciable rate gains, as compared to the maximum 
rate achievable by TDD operation. Intuitively, the condition 
states that a station’s amount of SIC should be at least as high 
as the loss incurred due to wireless propagation on the path 
to the intended receiver. 

Condition 5.4: ymm < and yi,i, < ' 

Proposition 5.5: If 7 mm < 1 +^’ ^ strictly 

concave and strictly increasing in Pm when Pi, is fixed. 
Similarly, if 755 < 


'TDD,m 


which is equivalent to ( 8 ) 




, r is strictly concave and strictly 
is fixed. Thus, when Condition 15.4 


increasing in Pi, when Pj 
holds, r is strictly biconcave and strictly increasing in Pm 
and in P 5 . Furthermore, when Condition |5.4| does not hold. 


r — Vr 


TDD 


< Ib/s/Hz. 


Proof: Fix P 5 . From the proof of Lemma mi we can 
express r as r = log(l + 7), where 7 is strictly increasing and 
concave in Pm whenever 

hmb^m ^ 

- 5 - -gm>0. (9) 

Nb + gbPb 

Multiplying both sides of (j^ by ^ and reordering terms: 


hmbPrr 


> 


gmPrr 


Nb + gbPb 


^ 7mm ^ 


5mb 


Ibb 


Whenever (|^, or equivalently, the inequality 7 mm < ’ 


holds, since 7 > 0 , > 0 , ^ 


< 0: 


dr 


1 


Q2rjr 


dy 

1+7 dPm 

97 N 2 


/ (31 y 
\dPm) 


> 0 , and, 

1 5^7 


< 0 , 


dPrn^ (1 + 7)2 \dP^) 1 + 7 dPj 

and therefore r is strictly increasing and strictly concave in 
Pm. Similarly, whenever 755 < , r is strictly increasing 

-L^r'Trn.m 

and strictly concave in P 5 when Pm is fixed. 

Now suppose that Condition |5.4| does not hold. Then, either 
Imm > 77 or 76 fc > 7^. Suppose that 7 ^^ > 7 ^. 
Then: 


r = 




1os(2-(i + 1(76™ + 


+ 766 

= Ib/s/H. + log (1 + i (7»„ + - 1 )). 

Since 1 (^^bm + ^ 76m}, it follows that 

r < Ib/s/Hz + which completes the proof for 7mm > 
. The proof for the case 755 > follows the same 

-L^r'Tbb -L^r'Tmm 

line of argument and is omitted for brevity. ■ 

B. Two Unidirectional Links 

Much of the analysis for use case (i) (Section |V-A| ) extends 
to use case (ii) (Fig. l||^b)| ), due to the similarity between the 
sum rate as a function of transmission power levels for these 
two use cases (see Eqs. 0 and 0). However, there are also 
important differences. First, the interfering signal at MS 2 in 
use case (ii), unlike the self-interfering signal at the MS in 
the bidirectional link case, is not known at the receiver, and 
therefore, cannot be cancelled (unless an additional channel 
is used, which we do not consider). Second, in use case 
(ii), the channel gains between MSs cannot take arbitrary 
values. This is because the channel gains typically conform 
to a path loss model of propagation, where the SNR depends 
on distances between MSs, which in turn need to satisfy the 
triangle inequality. The following two Lemmas are similar to 
Lemmas |53 and |5.3| We state them without proofs. 

Lemma 5.6: If there exists an FD sum rate that is higher than 
the maximum TDD rate, then the FD sum rate is maximized 
at Pmi = Pm for MS 1, and P 5 = P 5 for the BS. 

Lemma 5.7: FD extends the TDD capacity region by 
p • 100 % if and only if: 

.og(l + ^) .og(l + ^) 


log(l + 76^2 ) 


iog{i+7™:^) 


= l+p. ( 10 ) 


In a path loss model of propagation, the wireless channel 
gain between two stations is a function of the distance between 
the stations: huv = ( 7 ^) ’ where u^v ^ {6, mi, m 2 }, u v, 
T] is the path loss exponent, and P is a constant. Therefore, as 
distances dmib, dbm 2 ^ dmim 2 i^oed to satisfy the triangle 
inequality, SNRs ymib, lbm 2 and INR 7 mim 2 cannot take 
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Figure 7: TDD capacity region extension due to FD as a function of SNRs, where SNRs change due to path loss with exponent 
T], and distance between MS 1 and MS 2 is dmim 2 = p{dmih + Transmission power levels are set to maximum. In 

SNR regions where the triangle inequality of the distances is not satisfied, p is set to 0. 


arbitrary values. To evaluate rate gains in use case (ii), we 
consider path loss exponents r] G {2, 3,4}, since typical range 
for the path loss exponent is between 2 and 4 p0| . We assume 
fixed maximum power levels at the BS and the MS 1, equal 
noise levels N at the BS and the MS 2, and we vary SNRs 
and the INR as the function of distance, as follows: 

_ hmibPmi _ hmib ^ max _ f ^mib \ ^ max 

AT- Lmax imib I Jmin / imib’) 

'^mib 

_ ^mim2-Pmi _ ^m2m2 ^ max _ f ^mim2 \ ^ max 

7mim2 Ty- /,max imim 2 I ^min ) 7mim2’ 

'^77111712 '^mim2 

_ hi)rn 2 Pb _ hi)m 2 ^ max _ f ^bm 2 \ ^ max 
ibm 2 AT- ^max ibm 2 \ Jmin / ibm 2 ’ 

^bm2 ^%m2 

where is a reference distance at which juv = foi* 
u,v e {6,mi,m 2 }, X ^y. 

For the purpose of comparison, we will assume that = 
^mTb ~ ^m^m2 — which would Correspond to Pi, = Pm, 
and normalize all distances to d^^^. 

Capacity region extension as a function of SNRs is shown 
in Fig. 1^ for different values of the path loss exponent and 
17712 — pi^dmib H” di,m 2 )? P ^ }0.25,0.5,0.75,1}. For all 
combinations of SNRs at which the triangle inequality is not 
satisfied, we set the capacity region extension p to 0. 

Fig. [7] suggests that to achieve over 50% capacity region 
extension, the environment needs to be sufficiently lossy, i.e., 
with the path loss exponent 77 > 2. Moreover, to achieve high 
capacity region extension, the SNRs at the BS and at the 
MS 2 need to be low enough, meaning that the corresponding 
distances dmib and ^5^2 to be large, since the differences 
in the SNR shown in all the graphs are due to different 
distances (and consequently different path loss). 


VI. OFDM Bidirectional Links 


In this section, we focus on the rate maximization for use 
case (iii) (Fig. m- Recall that in this use case the FD receiver 
at the MS has a frequency-selective SIC profile (Fig. ^ ^ 


6.5 


Requiring two technical conditions (Conditions |6.2| and 
we derive an algorithm (Algorithm 1, MaximumRate) for 
the sum rate maximization. The algorithm is guaranteed to 
converge to a stationary point, which in practice is typically 
a global maximum. While the derived algorithm runs in 
polynomial time, its running time is high because it requires 
invoking a large number of biconvex programming methods. 
We therefore consider a high SINR approximation of the sum 
rate, and develop an efficient power allocation algorithm for 
the sum rate maximization. We also prove that in the high 
SINR regime it is always optimal to set the maximum SIC 
frequency in the middle of the used frequency band. 


A. Analysis of Sum Rate 


1) Dependence on Channel Power Levels: The analysis 
of the sum rate in terms of transmission power levels extends 
from the single-channel case (Section V-A[ ). In particular: 
Observation 6.1: If 


gmjk - cf 
Nm 


< 


dmb,k 


Nb + gbPb,k 


A db . 

and ^ < 


dbm,k 


P^b P^m P 9mPm,k{k c)^ 

( 11 ) 


hold, then the sum rate is biconcave in Pm,k and Pi,^k- 
This result is simple to show, since {k—c)‘^ term is independent 
of the transmission power levels, and Pi,^k and Pm,k only 
appear in one summation term (r/c). Therefore, we get the 
same form of partial derivatives in Pi,^k and Pm,k as in the 
case of a single channel (proof of Lemma [5T] ). Similar to the 
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case of a single channel, if condition is not satisfied, then 
the achievable rate improvement is low. 

The first inequality in (pTj) guarantees concavity in Pm,k 


when is fixed, while the second one guarantees concavity 
in Pi,^k when Pm,k is fixed. The condition (11) cannot be 


satisfied for any > 0, Pm,k ^ 0 the first inequality 

> However, since the 


cannot be satisfied if 


gmjk-c) 


_ _ _ Nm ^ Ni, 

role of condition ( pT] ) is to guarantee biconcavity in the power 
levels, we can replace this condition by either Pm,k = 0 
or Pi,^k = 0, which implies rate concavity in Pm,k^Pb,k- 
Specifically, to guarantee that the sum rate is biconcave in 
all Pm,k^Pb,k^ we require the following condition: 

gmjk-c)^ ^ hmb,k -f gmik-py 


Condition 6.2: (a) 


^b~\~gbPb,k 

Otherwise Pm,k = 0. and 


if 


Nrr. 


< 


(b)fr< 


hbrt 


_ if ^ ^ 


hbm,k . 


; otherwise Pb k = 


Nb - Nm+gmPm,kik-cy 

0 if Pm,k was not set to 0 by (a). 

Note that Condition 16.21 forces a channel k to be used in half¬ 
duplex (only one of Pm,k^ Pb,k is non-zero) whenever it is not 
possible to satisfy the sufficient condition ^TT\ for the sum rate 
biconcavity in Pm,/c, Pb,k for any Pm,k ^ 0 and P^^k ^ 0- 

2) Dependence on Maximum SIC Frequency: The 
following lemma shows that choosing optimal c for a given 
power allocation {Ph^k-, Pm,k} is hard in general, since the sum 
rate r as a function of c is neither convex nor concave, and 
can have Ct{K) local maxima. Proof is provided in Appendix. 

Lemma 6.3: The sum rate r is neither convex nor concave 
in c. All (local) maxima of r(c) lie in the interval (1,P). In 
general, the number of local maxima is Ll{K). 

Even though r(c) can have multiple maxima in c, if we 
restrict the analysis to the values of 7m6,/c and 7mm,/c that are 
relevant in practice, the selection of c, together with the power 
allocation, are tractable if the following inequalities hold: 


9m 


< 


hmb,k 


Nm Nb + 9bPb,k 


P}. 


( 12 ) 


Note that these inequalities are implied by Condition 6.2 for 
\k — c\ > 1, and that there can be at most 2 channels with 
\k — c\ < 1. For \k — c\ < 1, the corresponding inequality 
limits SI on channel k. The following lemma bounds the first 
partial derivative of r with respect to c. This bound will prove 
useful in maximizing r as a function of c and {Pb,k^ Pm,k} 
(Section rVEBT] ). 

Lemma 6.4: If inequalities ( p^ hold, then: 


dr 

dc 


< —iHK) + l + 2V3) Vce(l,i^). 


Similarly as for Condition 62^ since © cannot be satisfied 


for Pb > 0 when ^ > 


, we require the following: 


^mb,k '-C g-m ^ 
Nb+gbPb,k Nm 


Nm ^ Nb 

Condition 6.5: \/k G iT}: < 

, and Pm,/c = 0 otherwise. 

Proof of Lemma 6.4 can be found in the appendix. 


B. Parameter Selection Algorithms 

1) General SINK Regime: The pseudocode of the 
algorithm for maximizing the sum rate in the general SINR 
regime is provided in Algorithm 1 - MaximumRate. We 
claim the following: 


Lemma 6.6: Under Conditions |6.2| and |6.5| the sum rate 
maximization problem is biconvex. If biconvex programming 
subroutine in MaximumRate finds a global optimum for 
{Pb,k^Pm,k}^ then MaximumRate determines c and the 
power allocation {Pb,k^Pm,k} that maximize sum rate up to 
an absolute error e, for any e > 0. 

Algorithm 1 MAXiMUMRATE(e) 

Input: A, Pb, Pm 5 9b ? 9m ? Hm 5 A b 
1: Cl = 1, C 2 = K, Ac = _^(in(if)+i_,.2V3) 

2: c™’' = = 0, {PtXI = {-P™} = {0} 

3: for c = Cl, c < C 2 , c = c + Ac do 
4: Solve via biconvex programming: 

max r = where Vk is given by ^ 

s.t. Conditions 16.21 and 16.51 hold _ 

Y,k = l ^^^k < Pm, Pb,k < Pb 

Pb,k>0, Pm,k>0 ,V/CG A}. 

5: if r > then 

^max ^ ^max ^ 

V: {nX} = {PX} = {^-,4 

8 : return 


Note that without Condition 6.2 the biconvex programming 
subroutine in MaximumRate would not be guaranteed to 
converge to a stationary point (see, e.g., |[3T|). Moreover, 
since the sum rate is highly nonlinear in the parameter c 
(Lemma |6.3| ), c cannot be used as a variable in the biconvex 
programming routine (or a convex programming method). 
Nevertheless, as a result of Lemma |6.4| that bounds the first 
derivative of r with respect to c when condition |6. 5 1 is applied, 
we can restrict our attention to c’s from a discrete subset of 
the interval (IX)- 

Proof of Lemma 6.6' Consider the optimization problem 
in Step 1^ of the algorithm. Since Condition |6.2| is required by 
the constraints, the objective r is concave in Pb^k whenever 
Pm,k^ are fixed, and, similarly, concave in Pm,k whenever 
Pb^k^ are fixed. Therefore, r is biconcave in Pb,k^Pm,k- The 
feasible region of the problem from Step [^is determined by 


linear inequalities and Conditions |6.2| and 6.5 


Condition |6.2| is either an inequality or an equality for 
each Pm,k, Pb,k that (possibly rearranging the terms) is linear 
in Pm,k^Pb,k- Condition 6.5 is a linear inequality in Pm,k- 
Therefore, the feasible region in the problem of Step is a 
polyhedron and therefore convex. It follows immediately that 
this problem is biconvex. 

Suppose that the biconvex programming method from Step 
4 of MaximumRate finds a global optimum. Then the 
algorithm finds an optimal power allocation for each c from 
the set of _ 2 equally spaced points 

from the interval and chooses c and power allocation 

that provide maximum sum rate r. 

What remains to prove is that by choosing any alternative 
c and accompanying optimal power allocation the sum 

rate cannot be improved by more than an additive e. 

Recall from Lemma [63] that optimal c must lie in (IX)- 
Suppose that there exist c*, /d ^^‘^fi ^fi^^ ^ 

(IX), and r(c*, {P*^;PX}) + 

From the choice of points c in the algorithm, there must 
exist at least one point that the algorithm considers such 
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that I 

r(c 


c" — c 


< Ac = 


i^(ln(K) + l+2x/3)- 


From Lemma 


6.4 


: J) - r(c% J) < 

■ + 1 + 2\/3)^ = 


K 

E 

fc=i 


/log( W Wiog^ \y (13) 

^ ^ 1 “h 'ybb,k ' ^ 1 “h '~^m,m,.h ^ ' 


'^ram^k ^ 


While in the high SINR regime the dependence of sum rate 
on each power level Pb,kiPm,k foi* k G Ff} becomes 

concave (regardless of whether Condition [6^ holds or not), the 
dependence on the parameter c remains neither convex nor 
concave as long as we consider a general power allocation. 
Therefore, we cannot derive a closed form expression for c 
in terms of an arbitrary power allocation. However, as we 
show in Lemma |6.9[ when power allocation and the choice 
of parameter c are considered jointly, it is always optimal 
to place c in the middle of the interval c = 

The following proposition and lemma characterize the optimal 
power allocation for a given c. 

Proposition 6.7: Under high-SINR approximation and any 
power allocation {Pm,k} at the MS and any choice of c, it is 

'p 

always optimal to allocate BS power levels as Pb^k = 

Proof: Let Pb denote the total irradiated power by the BS. 
Write power levels on individual subchannels as = PkPb, 
where [3k > 0, V/c G {1,..., K}, and Ylk=i Pk = 1- Then the 
sum rate can be written as: 

n / hfYib kPm,k \ 


r = 


k=l 


+ 


log ( 


Nb + QbhPb' 

kbm,k[3kPb 


9m{k C)‘^Pm,k 


))■ 


First, observe that 
dr 

Qp^ ~ / ^ {/3fc>0} 

k=l 
K 

= E^{/3fc>o} ft 

k=l 


1 

n 


Nb 


QbPk \ 

+ gbf3kPb) 


-1 


9bl^k 


+Pb. 


dependence on each Pk^ and, moreover, is strictly concave 


in each /?fe, as ^ 




< 0, where 


i^(ln(if) + H-2x/3) 

since in any finite interval I any continuous and differentiable 
function f{x) cannot change by more than the length of the 
interval / times the maximum value of its first derivative f'{x) 
(a simple corollary of the Mean-Value Theorem). 

Since the algorithm finds an optimal power alloca¬ 
tion for each c, we have that r{c^, {P^P^ j^}) < 

r(c“, J) < r—. Therefore: r(c*, - 

^max ^ which is a contradiction. ■ 

2) High SINR Regime: A high SINR approximation of 
the sum rate is: 


the inequality follows from [3k < + [3k • Therefore, 

is maximized for [3k = p. ■ 

Lemma 6.8: Under high-SINR approximation and for a 
given, fixed, c the optimal power allocation at the MS satisfies 
Pm,k = Oik • Pm, where Ok >0, ^ak = 1, and for k ^ K: 


(i) “fe = (^ - 


-1 


(ii) ctk = 


Nm/RK+aK 


if = c, 


-j-4aKiNm+RKO!.K)Rk 


2Rk 


if k ^ c, 


where Rk = gm{k - crPm for k e {1, K}. 

Proof: Let Pm,k = Ok • Pm, where ak > 0, V/c, and 
Ylk=i^k = 1- The sum rate can then be written as: 

hmb,kOkPm \ 


r = 


K 

E 

k=l 


(log ( 


+ 


log ( 


p^b + 9bPb,k ' 

kbm,kPb,k 


:))■ 


^P'm 3~ Qmik 0)“^ Oik Pm 

Proving that at the optimal solution that maximizes r we 


necessarily ha ve P m, = Pj 
in Proposition 


6.7 


is analogous to the proof given 
^for Pi, = Pi,, and it is therefore omitted. 
As Y^k=i^k = 1, only K — 1 afs can be chosen 
independently, while the value of the remaining one is implied 
by their sum being equal to 1. Choose ok = 1 — ^k=i ^k, 
and observe that that K — c ^ 0 (and therefore Rk 0) is 
always true since, similar as in the proof of Lemma [63] at the 
optimum it must be c G (1, iT). 

If Rk = Qmfk — c)‘^Pm = 0, then the first and the second 
derivative of r with respect to ak are given as: 


dr 


1 1 doK 


dak 

Ok 

Ok dok 


1 

1 


= - 

+ 


Ok 

Ok Nn 

d^r 

1 

( 1 

dak^ 

“fe 

V'^lr (■ 
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( 1 


“fe 

Wk (■ 
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doK 


Nm + RrOk dOk 
Rk 


RkOk 
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{Nrr 


RrOk)^ 

1 


■OKf 
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It follows that r is concave in Ok and maximized for 

Ok = {ok~^ - {Nm/RK + Ok)~^) , 

where Rk = gm{K - cfPm- 

If Rk ^ 0, then the first and the second derivative are: 


(14) 


dr 

dok 

Q2r 


Rk 


- {Nb/{gbPk) P Pb) dokf 


Ok 
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P^i RkOk 


~ -2 ^ 

^k 


Rk 


Ok Nm 
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Rk 


RkOk ’ 
Rk^ 


{NmPRkOk)‘^ o\- {NmPRKOK)‘^ 


where is an indicator function. Since /5/c > 0 V/c G 
P} and Yl^=iPk = 1, it follows that there exists at 
least one strictly positive ^k- For each such [3k, — 

^ > 0, since Pb < + Pb. Therefore, ^ > 0, 


~ i^'m/Rk P Ok) ^ 

(f^K^ — {N^/Rk P Ok) < 0. 


— o 


which implies that it is optimal to choose Pb = Pb. 

Taking the first and the second partial derivative of r with 
respect to each [3k, it is simple to show that r has the same 


It follows that r is concave in Ok and maximized for: 

El = L_ ^ _L +__= 0. (15) 

dok Ok Ni -|- RkOk Ok Nm + RrOk 
After simplifying tig, we get: 

Ok{Nm P RkOk) = OK{Nm P RrOk). 


(16) 
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Solving the quadratic equation ( p^ for and using that > 
0 , it follows that r is maximized when satisfies 
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-Nm + Nm^ + AaK{Nm + RK(^K)Rk 
2Rk 
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lows that: maxr —r(c™^^, {P^ 


gm{k-cy 

- for Pm e [-Pto(1 - e'),Pm], it fol- 
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^Pm.e'. 


(17) 


where 7?/^ = gm{^ - cfPm^ ■ 

It is relatively simple to show (using similar approach as in 
the proof of Lemma |6.3| ) that under general power allocation 
r can have up to K local maxima with respect to c. However, 
if c is considered with respect to the optimal power allocation 
corresponding to 


c (Proposition |6.7| and Lemma | 6 . 8 | ), it is 
always optimal to place c in the middle of the interval ( 1 , 7f), 
as the following lemma states. 

Lemma 6.9: If {c, {Pi,^k, Pm,k}) maximizes the sum rate 




(b) K = {9,17,33} 


Figure 10: Evaluated (a) sum rate for K = 33, normalized to 


K, and (b) capacity region extension for K G {9,17, 33}, for 
the sum of the total transmission power levels at the MS and 
at the BS scaled so that it is the same as in the TDD case. 


under high SINR approximation, then c = . 

Even though this result may seem intuitive because the optimal 
power allocation is always symmetric around c (Proposition 
6.7 and Lemma | 6 . 8 | ), the proof does not follow directly from 
this property and requires many technical details. Eor this 
reason, the proof is deferred to the appendix. 

A simple corollary of Lemma |6.9| is that: 

Corollary 6.10: If (c*, maximizes r under 

high SINR approximation, then the power allocation 
is symmetric around and decreasing in \k — c\. 

Proof: The first part follows directly from 
The second part is proved in Lemma |6.9| ■ 

Lemma 6.11: A solution P™/?^}) that max¬ 

imizes r under high SINR approximation up to an absolute 
error e can be computed in O {K\og (^)) time. 

~ ~ ~ ‘ ‘ at the optimum P^” 


Setting: 


K 


• Pmc' = e e' = we yield the 


total running time of: O (Plog = O (Plog (^)). ■ 

We summarize the results from this section in Algorithm 2 
- HSINR-MaximumRate. 

Algorithm 2 HSINR-MAXiMUMRATE(e) 


Input. , Pi) , Pm 5 Qb ? Qm 5 Hm 5 Hb 

^max ^ ^K±1)/2 

{Pf:n = Pb/K. yke K} 

for ax G [0,1/A], via a binary search do 

Compute o/c for 1 < /c < A — 1 using ( pi) and ( p^ 
End binary search when G [1 — e/(A -h e), 

return 


Proof: Erom Proposition 

V/c G {1,...,A}. This can be computed in constant 
time, and requires 0 (A) time to assign the values to all the 
P 5 /c’s. Erom Le mm a 
Erom Lemma 6.8 


= akPm, where {ak} are positive 
coefficients given by (14), 0 and 'y OLk — 1* Rccal. that 
all the Qf/c’s are given in terms of ax, so we can find the 
allocation {ak} by performing a b inary search for ax until 


VII. Measurement-based Numerical Evaluation 

This section presents numerical evaluations for use case 
(iii). Numerical evaluations for use cases (i) and (ii) were 
already provided in Sections |V-A| and |V-B[ respectively. We 
focus on the impact of a frequency-selective SIC profile in a 
small form factor hardware at the MS (Eig. |^[d)| ), and evaluate 
achievable rate gains from ED. 


implies that ax < 


so it is sufficient to perform the binary search for ax C 
p] • a binary search requires 0 ( log {y^)) iterations, 
with each iteration requiring 0(A) time to compute [ak] and 
evaluate ^k, for the total time 0 (Alog 

The last part of the proof is to determine an appropriate 
e' so that > maxr — e, where the 

maximum is taken over all feasible points {c, {Pm,kj Pb,k})- 
Notice that we are only deviating from the optimal solution in 
that Ef=i = A* ■ Ef=i e [F^{1 - e'), Ai] instead 
of = Prn. Therefore, (c--, 

the optimal solution to the problem that is equivalent to the 
original problem, with maximum total power at the MS equal 
to Pm • Ylk=i ^k- Observe that: 

dr ^1 1 


Evaluation Setup. To determine the position of maxi¬ 
mum SIC and the power allocation R^k^} that maxi¬ 

mize the sum rate, we run an implementation of the Max- 
imumRate algorithm separately for measured ( p5j , p6| 
and Eig. |^{d)| ) and modeled (Eq. Q) SIC profiles of the MS 
ED receiver. Additionally, we determine {^m^k^R^k^} 
for the high SINR approximation of the sum rate using the 
HSINR-MaximumRate algorithm. We also compare the 
results to the case when the total transmission power is 
allocated equally among the frequency channels (we refer to 
this case as equal power allocation). 

Since the measurements were performed only for the ana¬ 
log part of the ED receiver, we assume additional 50dB 
of cancellation from the digital domainSimilar to [ p^ , 
we assume that when either station transmits at maximum 
total power that is equally allocated across channels (so that 
Rm,k = Rm/K^Rb,k = Rb/K)^ thc noise on each channel is 
llOdB below the transmitted power level. 

We consider a total bandwidth of: (i) 20MHz in the range 
2.13-2.15GHZ, (ii) lOMHz in the range 2.135-2.145GHz, 


Fig. ^ d) only shows isolation from the SIC in the analog domain. 
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Figure 8 : Power allocation over K = 33 channels (20MHz bandwidth) at the BS and MS for different values of average SNR 
(^avg) jYiQ higher the 7 ^^^, the more channels are used in full-duplex, and the closer the power allocation gets to the high 
SINR approximation one (computed by HSINR-MaximumRate). 



Figure 9: Evaluated |(a)| sum rate for K = 33, normalized to the number of channels K, and (b) - (d) capacity region extension 


for \(h)\ K = 33, = 17, and \(d)\ K = 9. The graphs suggest that higher average SNR ( 7 ^^^) and better cancellation 

(lower bandwidth - fewer frequency channels K) lead to higher rate gains. 


and (hi) 5MHz in the range 2.1375-2.1425GHz. We adopt 
the distance between the measurement points as the OFDM 
channel width 600kHz), so that there are K = 33, K = 17, 
and K = 9 channels, respectively, in the considered bands. For 
the SIC at the BS, we take = Ni, |l^ . 

We scale all the power variables so that Pm = P^ = 1. We 
consider flat frequency fading (so that hmb,k and h}jm,k are 
constant across channels k), and perform numerical evalua¬ 
tions for hmb,kPmlNb = hbm,kPb/Nm = * P, V/c, wtiere 

7 ^^g G {0,10,20,30,40,50} [dB]. 

We run MaximumRate for Ac = 0.01, which corresponds 
to an absolute error of up to e 0.2 for r. We evaluate the sum 

rate and the capacity region extension using the measurement 
data for the remaining SI and , {Pm^k ^ Pb^k^} returned 
by the algorithm. We assume that the amount of SIC around 
fc does not change as fc (and correspondingly c) is varied. To 
run the algorithm for c positioned at any point between two 
neighboring channels, we interpolate the measurement data. 
Results. Due to space constraints, we provide detailed results 
for the power allocation only for the 20MHz bandwidth (K = 
33) case, in Fig. For the lOMHz (P = 17) and 5MHz 


(K = 9) cases, we only provide the results for the capacity 
region extension, in Fig. 

Fig. shows the power allocations at the BS (Fig. |^[a)f[(d)| ) 
and at the MS (Fig. [(h)| ) computed by MaximumRate 
for both measured and modeled SI and for different values of 
average SNR 7 ^^^. Additionally, Figs, ^[d) and |^[h)| compare 
the power allocation computed by MaximumRate to the one 
computed by HSINR-MaximumRate. As Fig. suggests, 
when 7 ^^^ is too low, most channels are used as half-duplex 
- i.e., only one of the stations transmits on a channel. As 
yavg increases, the number of channels used as full-duplex 
increases: at 7 ^^^ = lOdB about seven channels are used as 
full-duplex, while for 7 ^^^ = 2OdB all but two channels are 
used as full-duplex, and when 7 ^^^ > 30dB, we reach the 
high SINR approximation for the FD power allocation. 

Fig. 1^ shows sum rate normalized to the number of 
channels for K = 33 (20MHz bandwidth) and |(b)f[(d)| capacity 
region extension for K = 33 (20MHz bandwidth), P = 17 
(lOMHz bandwidth), and P = 9 (5MHz bandwidth). As Fig.[^ 
suggests, the FD rate gains increase as 7 ^^^ increases and the 
SIC becomes better across the channels (i.e., as we consider 
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lower bandwidth - lower K). 

We observe in Fig. that there is a “jump” in the 

capacity region extension as 7 ^^^ increases from OdB to 5dB. 


This happens because at 7 ^^^ = OdB Conditions 6.2 and 6.5 


force all the power levels at the MS to zero, and we have the 
HD case where only the BS is transmitting. At 7 ^^^ = 5dB 
Conditions 16.21 and 16.51 become less restrictive and some of 
the channels are used as FD. At the same time, the total 
irradiated power (considering both MS and BS) is doubled 
compared to the case when 7 ^^^ = OdB (and to the TDD 
operation), so a large portion of the rate improvement comes 
from this increase in the total irradiated power. To isolate the 
rate gains caused by FD operation from those caused by the 
increase in the total irradiated power, we normalize the total 
irradiated power so that it is the same as in the TDD regime 
and compute the sum rate for K = 33 and the capacity region 
extension for K = {33,17,9}, as shown in Fig. 10 The 


results suggest that the rate gains that are solely due to FD 
operation increase smoothly with 7^^^ and the rate gains are 
almost indistinguishable for different power allocation policies 
(MaximumRate for measured and modeled SI, HSINR- 
MaximumRate, and equal power allocation). 

Since for the transmitted power of 1/K and c placed in 
the middle of the frequency band XINR at the first and the 
last channel is about 35 (?^ 15dB) for K = 33, about 8.5 
9dB) for K = 17, and about 2.5 (^ 4dB) for K = 9, 
our numerical results suggest, as expected (see e.g.. Figs. 
and Fig. [^, that to achieve high rate gains, 7 ^^^ needs to be 
sufficiently high. This is demonstrated by the results shown 
in Fig. and In particular, the rate gains obtained solely 
from FD operation are non-negligible when on most channels 
XINR > OdB. Moreover, simple power allocation policies, 
such as equal power allocation and high SINR approximation 
power allocation are near-optimal when the rate gains are non- 
negligible, as demonstrated by Fig. [T^ 


VIII. Conclusion and Future Work 

In this paper we considered three basic use cases of FD, 
including single- and multi-channel scenarios. In order to 
analyze the multi-channel scenario, we developed a new model 
that is grounded in realistic FD receiver implementations for 
small form factor devices. We characterized the rate gains in 
different scenarios and solved power allocation and frequency 
selection problems either analytically or algorithmically. Our 
numerical results demonstrate the gains from FD in scenarios 
and for receiver models that have not been studied before. 

This is one of the first steps towards understanding the 
benefits and the complexities associated with FD. Hence, 
there are still many open problems to consider. In particular, 
generalizing our results to the MIMO settings is of high 
relevance and interest. Additionally, SIC that has different 
impacts on different channels calls for the design of algorithms 
for OFDM networks with multiple access and MSs modeled 
as small form-factor devices. Moreover, we plan to develop 
scheduling algorithms that support the co-existence of half- 
and full-duplex users. While significant attention has been 
given to scheduling and resource allocation in half duplex 


OFDM networks (see, e.g., p2| and references therein), as 
demonstrated in this paper, the special characteristics of FD 
pose new challenges that have not been addressed. 
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Appendix 


Proof of Lemma 


6.3 


Since r = we will first 


observe partial derivatives of Vk with respect to c. 

Observe that in the expression ^ for Vk only ^rnm,k 
depends on c. Moreover, since ^mm,k = 
have that {k - c) = - 27 , 


gm{k-c)‘^Pm 
Nrr. 


we 


mm,k • 


Observe partial derivatives of with respect to c: 

dVk 2 QmPm,k 

~ \a2’ Nm ' 

k — c 

(1 ^mb,k ~\~ ^rnrn,k{cfj ^mm,k{6)^ 

_ 2 QmPm,k 

1 pj -J.J '^TTibjk 

In 2 Nm 

^mm,k{^Pj ^2 “h 'Tmb,k ^mb,k^ 


d‘^rk 

dc^ 


(1 'Tmbjk 


(18) 


(19) 


From dl 8 


^ equals zero for c = /c, it is positive for 
c <k and negative for c > k. Therefore, Vk is a has a unique 
maximum in c, with the maximum attained Sit k = c. Since 
this is true for every k G Ff}, it follows that for c < 1 

V/c G {1, ^ > 0 (with equality only for k = c), and 


therefore > 0. Similarly, ^ < 0 for c > iT. Therefore, all 
(local) maxima of r(c) must lie in the interval 

As 7mm,/c = ^ is symmetric around c= k. 

From (19), is negative for k — c = 0, and there exits 
a unique cq at which = 0 (this part can be shown by 
solving 7mm,/c (2 “h 7m6,/c 37mm,/c(f^)) (f “I“7m6,/c) — O? 


which is a quadratic equation in terms of {k—c)‘^ with a unique 
zero; see the proof of Lemma 6.4). For \k — c| > |A: - co|, 
is positive. This is trae, e.g., for -imm,k(c) > 1- 
Visually, each r/. as a function of c is a symmetric bell¬ 
shaped curve centered at k. Therefore, r can be seen as a sum 
of shifted and equally spaced symmetric bell-shaped curves. 
This sum, in general, can have linear in K number of local 
maxima. Examples with K local maxima can be constructed 
by choosing sufficiently large (sufficiently “narrow” 

bell-shaped curves). ■ 

of Lemma 6.4' Assume that jmm,k > 0 and 7m6,/c > 0 
\/k G {1, ...,Ff}, as otherwise |^| = 0 and can be ignored. 
Case 1. Assume first that c = /c* for some /c* G {1, ...,K}. 
Then, using ([^, ^ = 0, and for every k k"": 

QmPm 1^ c| 


dVk 


dc 
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m (1- T ^rr 
\k-c\ 


.,fc(c))(l -|“ ^mb,k^ 
^mb,k ^2 1 


In 2 1 -|- 3 m Pm t 

N rr 
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^-{k - c)2 1 +7m6,/c “ \n2\k- cY 

since k — c> 1. Observe that since c = /c* G {1,..., Ff}, every 
c — k is 3. positive integer. Therefore: 


dr 


drk 

< A 

k*-l K 

V ^ 1 V ^ 

dc 


k=i 

- In2 

7 —/c* ' \k — k*\ 

j=i ' fc=fc*+i ' ' 


- In2 L fc ^ 

k=l 


In 2 


where is the {K — 1)* harmonic number. Using the 

known inequality Hn < ln(n) +0.58+ ^ forn G N ||^ and 


assuming Ff > 4, we get: 1| < j^(ln(iT) +1). For K < A, 

by inspection: Ylk=i E ^ + 1- 

Case 2. Assume that c ^ iT}, and observe that for 

\k- 


> 1 : 1 % 


< 


< 


dc I — In 2 \k—c\ — In 2 [jk—c\ 

There can be at most two k's with 


'I'T- 


cl < 1. For such k, 


we bound |^| as follows. First, observe from (18) and 0 


d\k—c\ 

if for some cq: 


that^||;^| = - 


d Vk 
dc^ 


From 


0, 


= 0 if and only 


,fc(co) (2 + ^mb,k + ‘3^mrn,k{eff) ^mb,kY — 0 

(2 + ^mb,k) + \/(2 + 'ymb,k)‘^ + 12(1 + ^mbf) 
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^^7'nim,fc (co) — 

Note we have used that ^mm,k > 0 to get a unique solution 

for ^mm,k' Since 7nTm,/c(co) — ’ if Cq) . 

(fc - co)^ = + -7inTn,fe(co) 


> 

From condition 

gives \k - Col 


9mPm,k 

Nm 2 • (2 + 'Jmbjk^ 




> 


Cm,k O 

we have that 


Nrr 


6.5 vv^ C..C.C p 

9 m m, k 

> It is clear from 


Nm 'Jmbjk 
QmPm,k 3 
• lmb,k > h which 


and ^mm,k — 

^^{k — c)‘^ that is negative for \k — c\ < \k — co\ and 


Nm dc^ 

positive for \k — c\ > \k — co\. Since 


d 

d\k—c\ 


follows directly that I I is maximized at c 
for \k 


I ^Lk I _ 
I dc I 


d rk 
dc^ ’ 


it 


< 1, we have that | ^ | 


< 


Cq. Therefore, 

< +73- 


In 2 \k—co I ^ In 2 

Combining the results for \k — c\ > 1 and \k — c\ < 1: 


< + 1^1 < J_[|_ + 1 

+ lacl ln2ll ^ li- 

K-1 . 

{T.I 


+ 


L 

k=[c\ -1-1 


\k-c\ 


+ 2 


ln2 ’ 


+ 2 


< + 1 + 2a/3). 
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C I 



(a) c — 5 < I 


(b) c - 5 > I 


Figure 11: Pairing of points that are left and right from c for 
(5,5.5) andl^cG (5.5,6). 


of Lemma 


6.9 


From Lemma 


6.8 


^,/e = ^kPm. where 


Pm — Pmf 3.nd recalling that R]^ — QmPmiJ^ P) • 

Oi KK QmP m(-^ c) ) ./> 7 
• OLk = -^- N - - lfk = C\ 


_ -Nm + yJ Nm‘^+A(XK{Nm+CXKgmPm{K-cy)gmPm{.k-c)‘^ 
— 2gmPm{k-cy 

\i k C\ 


and ax > 0 is chosen so that Ylk=i <^/c = 1- To simplify the 
notation, we will let ^rnm = , and write ak as: 

foK • (l + aKjmmiK - c)^), if /c = C, 

“ I -l + \/l+4aj^(l + ax7mm(^-c)^)7mm(fe-c)^ if A; C ' 

( 20 ) 

Notice that for c = 1 + / • / G {l,2,...,2iT — 3}, the 

power allocation is symmetric around c, that is:Q^L|j = Q^[|], 

The first partial derivative of r with respect to c is: 

1 


dr _ 6 I 1 

gc~^ nr [ 


k=l 

K 


iV,n(l+ 2^afc(fc-c)2) 


k=l 

K 


^ dc \1 +7mmOfc(A; - c)2 

27 

mm akjk - c) 


= E 


■ 7 ^ 1 T ^mmOLkik c)^ 


( 21 ) 


Fix any I G {0,l,...,iT — 2} (on the left half of the interval 
[1, K]) and let c G (1 + / • |, 1 + (/ + 1) • ^). We make the 
following three claims: 

(Kl) Each point i G {1,2,..., [cj} (left from c) can be paired 
to a point j G {[c], [c] +1,..., Ff} such that all the pairs 
are mutually disjoint and for each pair (i, j) we have that 
c -i < j - c. 


Proof of (Kl): To construct the pairing, observe that, by the 
choice of c, c is between two consecutive integer points and 
is strictly closer to one of them. If it is closer to the left point, 
then the pairing is ([cJ, [c]), ([cJ - 1 , \c] +1),..., ( 1 , 2 [cJ). If 
c is closer to the right point, then the pairing is ([cJ, \c] + 1 ), 
([cJ — 1 , \c] + 2 ),..., (l, 2 [cj + 1 ). Such pairings must exist 
as c < The pairings for K = 12 and cases: c G (5, 5.5) 
and c G (5.5,6) are illustrated in Fig. 11 Q.E.D. 


(K2) In the optimal power allocation that corresponds to a 
given c and for any i^j G {1, ...,iT}, if \i — c\ 
then ai > aj.ln other words, the smaller the distance 
between k G {1,..., Ff} and c, the larger the ak- 


Proof of (K2): The proof has two parts. First, assume that 
|i — c| = 0 and observe aj for \j — c| > 0. From (20): 

oti = axil + OLKlmm^K - c)^), and 


1 T c) ^)'yrnm (j c) ^ 

27mm (j-c) 2 

1 T s/l T 4oi7mm(j c)^ 

27 mm (j c)^ 

Using simple algebraic transformations: 

OLj ■<C OLi 

1 T 1 T 4Q;j7mm(j c)^ 

27mm (j-c) 2 ^ 

^^s/l T 4Q;i7mm(j c)2 < 1 -f 2Q;^7mm(j c) , 


Observe that given the optimal power allocation ( [^ : 

• If c = then (from (20)) ai = ax, OL2 = OLK-i^- 


= a^K±iy and it follows that ^ = 0 . 


If c 
then, 

>LcJ _2ri 

mm 


— 1 + / • 

as {ttfe} is 


ai(i-c) 


E7 


X^LCJ _ 

1 I77mm ('i c) 

E K 2^rnm<^k jk — c) 

/c=2LcJ+l l-k'ymmOikik-cP 

E K 2jrnm0^kik — c) 

/c=2LcJ+l l+^^^afc(/c-c)2 

If c = ! + ;• i, for / e {K,...,K 
symmetric around c: 

j jj c) 


dr 
dc 

for I G {0,1,..., K - 
symmetric around c: 

2LcJ 

mm o^jU-c) 


[cj+l l-\-JrnmOij{j-cP 


2 }, 


+ 


> 0 . 


2 }, then, as {ak} is 


_ 

Z^j=2c—K 1 + 7 , 
^-^2c—1 
2^i=l 


mm 

mm 


ocfj-c) 

cxi{i—c) 


-c) I 

-c)2 ^Z^/c = 


1+7 mm 
2 + 


ai(i—c)2 ^ 

akjk-c) 


[cj+l 1+7 

mm akik-cP 


< 0 . 


l+7mm('^ c)^ 

In other words, if we restrict our attention only to those [ak] 
that determine the optimal power allocation, then considering 
c’s from the set 1 + / • |, where I G {0,1,..., 2K — 2}, we 
get that the first derivative of r with respect to c is positive 
for c < equal to zero for c = and negative for 

c > To conclude that at the global maximum for r we 

have c = by considering c G (1, Ff) it remains to show 
that for c G (l+Z-^, l + (/H-l)'|), where I G {0,1,..., 2K—2}, 
we have that > 0 if / < - 2 and < 0 if / > - 1 . 


we get that aj < a^ by squaring both sides of the last term, 
as \j - c| > 0 implies {2aiyrnm{j - c)^)^ > 0 . 

Second, assuming that \k — c\ > 0 and taking the first 
derivative of ak with respect to (k — c)‘^, we show that ak 
decreases as (k — c)‘^ (and consequently \k — c\) increases. Let 
A = (k — c)^. Then, as: 

dOLk d ( 1 "y/1 + ^OLxil T OLxymmi^K C)‘^)ymmA 

dA dA y27mm^ 27mm ^ 

_ 1 _ 1 + 2ok(1 + axy mm {K - c)^)7mmA 

27 mm 27 mmA 2 y^l + 4 ok(1 + OiKymm(K - c)2)7 mm A 


it follows that ^ < 0 , since _ 

1 + 2aK{l + aK-fmmiK - c) 2)7 • Q-E.D. 

(K3) As \k — c\ increases, decreases. 

I OC I i-t-7mmQ^fc CJ 

Proof of (K3): Observe that: 


d dvi^k 
dak dc 


27mm I ^ c\ 

(1 + 7mmOfc(/c - C)2)2 


> 0. 


We had from 


(K2) 


that 


dcx k 

d\k—c\ 


< 0 , and therefore: 


d 

dvk 

d 

dru 

dak 

d\k — c\ 

dc 

dak 

dc 

d\k — c\ 


< 0, Q.E.D. 
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Using (|21|, we can write ^ as: 


K 


dTj _ 
dc ^ 


c) 

1 ~\~ 'ymm0^k(^k c)^ 

LcJ 


_ ‘^Jmm 0 ^i{i — c) ^ 

" ^ l+7mm«i(i-c)2 


^'ymmO^ji^j c) 

1 H“ ^mmOLj (j c)^ 


If c G [1, then, from (Kl) 


each term i in the left 
summation can be paired to a term j in the right summation, 
such that all the pairs are disjoint and for each pair (i, j): 
H — c\ < \j — c\. From (K3) for each such pair (i, j): 


1'^ c| 


mm^j (^3 

left summation are negative, and all the terms in the right 
summation are positive, it follows that: 

2 ' 7 mmOi(z c) c) 


I J C\ 


As all the terms in the 


dr 


uf 

dc 1 


+ 7rr 


i{i - c)- 


+ 


E 


LcJ 


=-5:t 


27rr 


1 H“ ^mn 
27rr 


7(i -c)2 


1 “1“ c)‘ 


j=[cj+l 
K 

j=^+l 1 + lmmaj{j - C)2 


^aj|i - c| 


> 0. 


Proving that < 0 for c G is symmetrical to the 

proof that |^ > 0 for c G [1, As |^ = 0 for c = 

at the globally maximum r we have that c = . ■ 













